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SUMMARY 

The dynamic behavior of the wave journal bearing was determined by running a three-wave bearing with 
an eccentrically mounted shaft. A transient analysis was developed and used to predict numerical data for the 
experimental cases. The three-wave journal bearing ran stably under dynamic loads with orbits well inside the 
bearing clearance. The orbits were almost circular and nearly free of the influence of, but dynamically 
dependent on, bearing wave shcpe. 

Experimental observations lor both the absolute bearing-housing-center orbits and the relative bearing- 
housing-center-to-shaft-center oibits agreed well with the predictions. Moreover, the subsynchronous whirl 
motion generated by the fluid film was found experimentally and predicted theoretically for certain speeds. 


SYMBOLS 

B damping coefficient between bearing housing and its support, N s/m 

B ijii = rrJ=rn film damping coefficients, N-s/m 
C bearing radial clearance, m 

F,., F, fluid film force ccmponents along and perpendicular to line of centers - O (fig. 1 ), N 

F i(Xi= rn steady-state fluid film force components along and perpendicular to line of centers - O, N 

F;, F ^ fluid film force components along and perpendicular to line O 0 - O (fig. 1), N 

h fluid film thickne >s, m 

K stiffness coefficient between bearing housing and its support, N/m 

K ljii= rrJ = rJ) film stiffness coefficients, N/m 
M bearing-housing total mass, kg 

O bearing-housing center (fig. 1) 

O 0 fixed center of roiation (fig. 1) 

O, shaft center (fig. i ) 

p fluid film pressure. Pa 

R shaft radius (bearng normal radius), m 

r coordinate along ine of centers 

S space (displacement) between shaft center Oj and bearing-housing center O, m 

S iU = rn S components along and perpendicular to line of centers 0 } - O, m 

t time, s; coordinate perpendicular to line of centers 

V velocity between shaft center Q ] and bearing-housing center O, m/s 


NASA/TM— 2002-211079 


1 



V iUs ,. n V components along and perpendicular to line of centers O, - O, m/s 

V n difference between shaft surface and bearing surface speed projected on perpendicular direc- 

tion to shaft surface, m/s 

Vy component of shaft surface speed along its circumference, m/s 

- axial coordinate parallel to shaft axis 

£ bearing eccentricity relative to shaft (relative movement) (fig. 1), 0 { - O, m, also bearing 

line of centers 

e/C bearing sleeve-to-shaft eccentricity ratio 

£ 0 bearing run-out (absolute movement) (fig. 1), O 0 - O, m 

£o/C bearing-housing absolute eccentricity ratio, £ 0 = O 0 - O 

£, T| axes along and perpendicular to direction O 0 - O (fig. 1) 

0 angular coordinate along shaft circumference, rad 

p fluid film dynamic viscosity, N-s/nr 

p fixed shaft run-out (fig. 1), O 0 - Oj, m 


<p angle between line of centers Oj - O and O 0 - O (fig. 1 ), rad 

\|/ rotation angle of O 0 - O around O 0 (fig. 1), rad 

Q rotation angle of O 0 - O, around O 0 (fig. 1), cor, rad 

co angular rotation speed (fig. 1), rad/s 


INTRODUCTION 

The wave bearing concept has been under development since 1992. Thus, the steady-state and dynamic 
performance under fixed side load (refs. 1 and 2) and the influence of both the number of waves and the ratio 
of wave amplitude to radial clearance (refs. 3 and 4) have been analyzed. Moreover, the steady-state charac- 
teristics of the wave journal bearing and its dynamic stability have been experimentally measured. Good 
agreement was found between the experimental data and theoretical predictions (refs. 5 to 8). In addition, 
the experimental work revealed good dynamic behavior of the wave bearing when subsynchronous whirl 
motion occurred. The wave bearing performed well, keeping the orbit of the subsynchronous motion inside 
the bearing clearance (refs. 7 and 8). Consequently, the wave bearing should perform well under the dynamic 
loading conditions that often occur in most rotating machinery. Any rotor can be subject to a dynamic load 
caused by an unbalance, or a run-out, of the shaft. This dynamic load is a rotating load that has a rotational 
speed equal to the rotor speed. Such a load can be simulated by running the bearing with a shaft that has a 
fixed run-out. Therefore, a transient analysis was performed to predict bearing behavior under a rotating load. 
Then an experiment was conducted to record the orbits of the bearing-housing center when the shaft has a 
known fixed run-out. 


ANALYSIS 

Bearing-housing-center movement can be studied by using the motion equation of the center along and 
perpendicular to the radial direction O 0 - O (axes £ and r| in fig. 1): 
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Also, figure 1 shows that the eccentricity e = O, - O (where O, - O is the line joining the shaft center O, and 
the bearing-housing center O) and that the shaft run-out p = O 0 — Oj. Assuming that the motion starts from the 
downward vertical where the shaft and the bearing are concentric (£ = 0), when the shaft rotates around O 0 
with the angular speed co, p makes the angle SI and drives the bearing so that £ 0 makes the angle vj/. 

The governing equations (1) are two scalar, coupled, nonlinear ordinary differential equations. These equa- 
tions are integrated simultaneously by using a fourth-order Runge-Kutta method for known values of M, C, T-, 
F n , K , and B and initial values of £ 0 , V|/, dzjdt, and d\\)/dt (ref. 9) The fluid film forces applied to the bearing 
surface are 


F- = F r cos tp+F, sin cp 

( 2 ) 

/-jj = F r sin cp - F, cos cp 

The projections of fluid filr l force along and perpendicular to the line of centers are 


F, = F r<] + K rr S r + K n S, + B lr V r + B n V r 
F r = F, 0 + K lr S r + K„S, + B„V r + B„V, 


<3> 


The bearing steady-state fo re and dynamic stiffness and damping coefficients can be computed by inte- 
grating the Reynolds pressure equation at each time step location of the shaft with respect to the bearing. This 
equation, assuming the gas will expand isothermally, is 
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(4) 


The Reynolds equation (4) can be integrated by using its complex form and a small perturbation technique. 
This procedure is described, foi instance, in reference 10. 

The solution procedure can start with an input data set (bearing length, diameter, radial clearance, shaft 
turning speed, shaft run-out, and the time step). In addition, a set of starting values at time = 0 are required: 


= P* 


dt 


V = 0, 



(5) 


Then, at each time step, where £,„ \p, and Q (Q = (M ) are known, the O 0 O,O triangle (fig. 1) is known, and 
all geometrical parameters as v'ell as displacements and velocities can be calculated. Therefore, the Reynolds 
equation (4) can be integrated over the fluid film. Then, all parameters of the motion equation ( 1 ) are known 
as well as the starting values fcr the next time step (E 0 , \p, and their time derivatives £2). The procedure is 
repeated until the orbits are completed. 
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APPARATUS 


The wave bearing rig described in references 5 to 8 was used to perform the experimental work. The axis 
of the spindle that drives this rig is vertical, and the experimental bearing housing is mounted on the rig table 
and supported by two pressurized thrust plates. This configuration keeps the bearing housing stiff in the axial 
and angular directions but allows it to move freely in the radial direction. The experimental shaft is an exten- 
sion of the rig spindle shaft. It is mounted into the tapered end of the spindle shaft with a fixed run-out (for this 
experiment 1 1 ±0. 1 pm). A cross section by a horizontal plane of the experimental bearing is shown in figure 2. 
The fixed rotation center for the system is O 0 . The centers of the shaft and the bearing housing are O, and O, 
respectively. The shaft run-out O 0 - O, is fixed. 

The goal of this work was to record and predict the absolute and relative orbits of the bearing-housing 
center O. The motion of the center O can be observed like an absolute motion for instance with regard to 
the center of rotation O 0 or like a relative motion with regard to the center of the shaft Oj. Figure 3 shows the 
experimental bearing setup. Two sets of light-beam proximity probes were used. Two probes were located at 
90° in the bottom side of the bearing housing and “looking" at the shaft. These probes detected the orbit of 
the bearing-housing center relative to the shaft center (O - O,). The second set of two probes were located 
also at 90° but held by supports fixed on the rig table and “looking" at the bearing housing. These latter probes 
detected the absolute orbit of the bearing-housing center (O - O 0 ). A polished circumferential strip was made 
on the outside bearing-housing surface to avoid asperity noise from its roughness. The light-beam probes were 
calibrated by using the known fixed run-out of the shaft. The displacement of the shaft was measured with a 
precision of 0.1 pm. The theoretical predictions of the orbits were made through a transient analysis of the 
bearing-housing-center motion. 


RESULTS AND DISCUSSION 

The experimental bearing was 5 1 ±0.01 mm in diameter, 58±0.01 mm in length, 20±1 pm in radial clear- 
ance, and 2.2±0.01 kg in mass. The bearing had three waves with a 0.5+0.07 ratio of wave amplitude to radial 
clearance. The shaft was set with an 1 1 ±0. 1 pm run-out. The damping, B in eqs. (1), in the bearing-housing 
support and connection system was found to be 0.05 N-s/m. The stiffness, K in eqs. (1), had little influence on 
the bearing orbits and was approximately zero. The top proximity probes (fig. 3) produced 500 mV for 5.78- and 
4.78-pm displacements in the horizontal and vertical directions, respectively, and the bottom probes produced 
500 mV for 6.11- and 6.90-pm displacements in the horizontal and vertical directions, respectively. (Horizontal 
and vertical directions refer to the directions on the oscilloscope photographs shown on the right sides of fig- 
ures 4 and 5, 90° apart in the physical plane. 

The test rig was run at four speeds up to 5540 rpm. Below 3100 rpm both the observed and predicted orbits 
of the bearing-housing center showed that a subsynchronous whirl motion took place inside the bearing clear- 
ance. Figure 4 shows both the predicted and observed orbits for relative and absolute motion of the bearing- 
housing center when the shaft rotated at 2156 rpm. When the speed increased above 3100 rpm, the motion 
stabilized, as shown in figure 5 for a shaft speed of 5539 rpm. 

Both the absolute and relative observed orbits of the bearing-housing center are shown as oscilloscope 
photographs on the right sides of figures 4 and 5. On the left sides of these figures the computed orbits are 
presented with a time step of 0.000001 s (10 ps) and for 30 000 steps. The experimental orbits appeared as 
ellipses rather than circles because of the difference in the probe sensitivity in the horizontal and vertical 
directions mentioned above. Both experimental orbits in figure 4 have a specific pattern caused by subsyn- 
chronous whirl motion. The transient analysis also revealed this pattern. Both the experimental and theoretical 
absolute orbits (fig. 4(a)) were within a radius of 5 to 12 pm. Both the experimental and theoretical relative 
orbits (fig. 4(b)) were within a radius of approximately 5 pm. 

The bearing stability increased as the running speed of the rig increased. Figure 5 shows the results for 
5539 rpm. The experimental orbits were perfectly stable. The shaft run-out made large absolute orbits of the 
bearing housing (right side of fig. 5(a)). However, the radius of the relative orbits was approximately 2.5 pm 
(right side of fig. 5(b)) despite the 11±0.1 pm shaft run-out (i.e., the bearing followed the shaft very well). The 
predicted orbits, shown on the left side of figure 5, matched very well with the observed orbits. The theory also 
showed that the bearing would run stably. After a couple of rotations from the starting point the orbits were 
stable, keeping almost the same path. 
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The relative orbits of the bearing housing increased but the absolute orbits decreased as the speed 
increased. This effect showed the influence of both external damping and bearing inertia on the magnitude 
of orbit radius. In addition, the rearing actually ran more and more stably as speed increased, and the theory 
showed that the number of rotations before the bearing reached a stable orbit would diminish as speed 
increased. 

All runs showed only a smt 11 influence of the bearing wave shape on the orbit shape despite the experi- 
mental bearing’s large wave amplitude ratio, 0.5±0.07. This result confirmed that a wave bearing with few 
waves, such as three, w'orked well under dynamic loads. The bearing behaved in such a way as to average the 
influence of the waves. 

Two types of shaft-centerec motion can be defined with respect to the center of the bearing: (i) stable 
unbalance or run-out orbits (e.g , fig. 5(b)), where the center of unbalance rotates at shaft frequency; and 
(ii) fluid-film-induced unstable whirl orbits (e.g., fig. 4(b)) that are superimposed over the stable unbalance or 
run-out orbits at a specific frequency different from the rotation frequency. The unbalance motion (i) is seen in 
each graph, but the unstable whirl (ii) occurs only at specific rotational speeds. 


CONCLUSIONS 

The dynamic behavior of the wave journal bearing was determined by running a three-wave bearing with 
an eccentrically mounted shaft. The following conclusions were reached: 

1. A dynamically loaded three-w'ave journal bearing can run stably, averaging its behavior when the 
wave exposure to the load is changing. The orbit radius of the relative motion between the shaft and 
the sleeve is smaller than the bearing clearance, and the motion is contained within the bearing 
clearance. The orbits aie almost circular and nearly free of, but dynamically dependent on, the influ- 
ence of bearing wave shape. 

2. Good agreement between experimentally observed and theoretically predicted orbits was found at all 
tested speeds for both relative and absolute motions. 

3. The subsynchronous whirl motion influences the bearing-housing-center orbits if the bearing speeds 
are in the region where the bearing itself is susceptible to subsynchronous whirl instability. When the 
bearing runs under such circumstances, the orbits show a specific pattern. This pattern was observed 
experimentally and was also confirmed theoretically by the transient analysis. 
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Figure 1 . — Geometry of three-wave 
journal bearing. Rotation center, Oq; 
shaft center, Oi ; and bearing-housing 
center, O. 



Figure 2. — Cross section of three-wave 
journal bearing by horizontal plane. 



Figure 3. — Experimental bearing setup. 
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Figure 5.— Predicted and experimentally observed orbits of three-wave journal bearing at 5539-rpm shaft 
rotating speed, (a) Absolute bearing-housing-center orbits, (b) Relative bearing-housing-center-to-shaft 
center orbits. 
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